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ABSTRACT 

Context. 

Aims. We present a new grid of massive population III (Pop III) star models including the effects of rotation on the stellar structure 
and chemical mixing, and magnetic torques for the transport of angular momentum. This grid covers the range of mass from 10 to 
1000 M Q , and rotational velocity from zero to 100% of the critical rotation on the zero-age main sequence. Based on the grid, we also 
present a phase diagram for the expected final fates (i.e, core-collapse supernova, gamma-ray bursts and pair-instability supernovae 
of diverse types) of rotating massive Pop III stars. 

Methods. The model grid has been calculated with a stellar evolution code. We adopted the recent calibration made with the VLT- 
FLAMES data for the overshooting parameter and the chemical mixing efficiency due to rotation. The Spruit-Tayler dynamo was 
assumed for magnetic torques. 

Results. Our non-rotating models become redder than the previous models in the literature because of the larger overshooting pa- 
rameter adopted in this study. In particular, convective dredge-up of the helium core material into the hydrogen envelope is observed 
in our non-rotating very massive star models (> 200 M Q ), which is potentially important for the chemical yields. On the other hand, 
the stars become bluer and more luminous with a higher rotational velocity. With the Spruit-Tayler dynamo, our models with a suffi- 
ciently high initial rotational velocity can reach the critical rotation earlier and lose more mass as a result, compared to the previous 
models without magnetic fields. The most dramatic effect of rotation is found with the so-called chemically homogeneous evolution 
(CHE), which is observed for a limited mass and rotational velocity range. CHE has several important consequences: 1) Both primary 
nitrogen and ionizing photons are abundantly produced. 2) Conditions for gamma-ray burst progenitors are fulfilled for an initial mass 
range of 13 - 84 M . 3) Pair instability supernovae of type Ibc are expected for 84 - 190 M e . 4) Both a pulsational pair instability 
supernova and a gamma-ray burst (GRB) may occur from the same progenitor of ~ 56 - 84 M Q , which might significantly influence 
the consequent GRB afterglow. We find that CHE does not occur for very massive stars (> 190 M Q ), in which case the hydrogen 
envelope expands to the red-supergiant phase and the final angular momentum is too low to generate any explosive event powered by 
rotation. 

Key words. Stars: evolution - Stars: rotation - Stars: Population III - supernovae: general - Gamma-ray bursts:general 



1. Introduction 

The first generations of stars (or population III stars) in the 
early Universe are supposed to play a key role in shaping the 
early Universe. It is believed that the first stars were intrinsi- 
cally massive, mainly because of the lack of efficient coolants in 
the primordial gas (e.g.|Bromm et al. ||1999 |Abel et al. |2000 



Nakamura & Umemura 200 1J. This means that they should be 
one of the most important sources of ultra-violet photons for 
reionization at high redshift. The chemical evolution of the early 
Universe is critically determined by the first stars because most 
of the elements heavier than helium are uniquely provided by 
them via supernova explosions and possibly stellar winds. The 
chemical abundance patterns observed in extremely metal-poor 
stars or in damped Lyman-a systems at high redshift might im- 



print the nucleosyntheses in the first stars (e.g. Depagne et al. 


2002 


Christlieb et al. ||2002 Cayrel et al. 2004| Erni et al. 


2006| 


Frebel & Norris |201 1 |Cooke, R. et al. |201 1 | |Kobayashi 


et al. 


2011b. Some of the first stars may produce very bright 



events such as supernovae (SNe) or gamma-ray bursts (GRBs) 
which are potentially observable. The recent discoveries of some 
GRBs at very high redshift ( Kawai et al. ||2006 Greiner et al. 
2009] ISalvaterra et al. ||2"0091 |Tanvir et al. ||20091 |Cucchiara et 
al. |201 l| l indeed imply that such luminous events from the first 



stars may serve as indispensable tools for the probe of the early 
Universe with the next generation of telescopes. 



Many authors already addressed these questions using stel- 
lar evolution models (Woosley & Weaver~|[ 1995; Marigo et al. 



200"T]|Heger & Woosley ||2002| |Marigo et al. ||2003| |Chieffi & 
Limoiigi |2004||Umeda & Nomoto |2003||2005||Ekstrom etaT 



2008 ; He ger & Woosley |2010[ ). One of the most intriguing ques- 
tions here is how rotation affects the evolution and the final fates 
of the first stars. Previous studies on massive star evolution indi- 
cate that both the rotationally induced chemical mixing and the 
centrifugal force can significantly alter the mass loss rate, the 
stellar structure, and the consequent evolution and nucleosyn- 
thesis ( |Maeder & Meynet ||2000b| |Heger et al. ||2000] >. This ef- 
fect should become even more important for metal-poor massive 
stars (e .g. |Meynet & Maeder |2002} |Yoon, et al. |2006| |Hirschi| 
1 |2007| l, for which angular momentum loss through radiation 



driven stellar winds is generally less efficient (e.g., IVink et al. 



2001 1 Kudritzki 2002). In addition, a recent numerical study 



on star formation in the early Universe by Stac y et al. | ( |201 1) 
shows that the rotational velocity of the first stars can be close to 
the break-up limit, suggesting that the effects of rotation on the 
evolution of Pop III stars should not be ignored. 
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There exist only a few studies on the role of rotation for the 
evolution of the first stars. Marigo et al. ( 2003| l calculated rigidly 
rotating massive metal free stars (120 - 1000 M Q ), finding that 
the centrifugally driven mass loss may be significant when the 
rotational velocity at the equatorial surface reaches the critical 
limit. Ekstrom et al. ( 2008| > considered the internal transport 
of angular momentum due to hydrodynamic processes and rota- 
tionally induced chemical mixing in their evolutionary models 
for masses of 9 - 200 M Q at zero metallicity. They showed the 
importance of the centrifugally driven mass loss, as well as the 
interesting consequences of chemical mixing in the nucleosyn- 
thetic yields. 

In the present paper, we present a new grid of evolutionary 
sequences of metal-free massive stars including the transport 
of angular momentum and chemical species, for initial masses 
ranging from 10 M G to 1000 M . While only a fixed value for the 
rotational velocity was assumed in Marigo et al. and Ekstrom et 
al. we considered several different rotational velocities for each 
mass. As shown below, this enables predicting the final fates of 
the first stars (core-collapse, GRB, pair-instability SN, etc.) in 
the parameter space spanned by the initial mass and the initial 
velocity (Sect. IS). The present study also differs from the pre- 
vious work because we included magnetic torques according to 
the Spruit-Tayler dynamo (Spruit 2002} . The adopted efficiency 
of the angular momentum transport in our work is therefore in- 
between that of Marigo et al., who enforced rigid rotation, and 
that of Ekstrom et al., who considered only hydrodynamic trans- 
port processes (see Sect. [2] for a more detailed discussion on 
this). 

The paper is organized as follows. We describe the numer- 
ical method, physical assumptions, and the newly constructed 
model grid in Sect. [2] The evolution of basic physical properties 
for different initial masses and rotational velocities and its con- 
sequence on the ionizing fluxes from rotating Pop III stars are 
discussed in Sect. [3] The effect of chemical mixing due to con- 
vection and rotationally induced hydrodynamic instabilities on 
the nucleosynthesis is investigated in Sect. |4] The conditions for 
the so-called chemically homogeneous evolution are discussed 
in Sect. [5] In Sect. [6j we explain how the final structure of the 
first stars would be affected by rotation, and discuss the possible 
diversity of the explosions of the first stars. Finally we summa- 
rize our results and discussions in Sect. [7] 



2. Numerical methods 

2.1. Basics on the stellar evolution code 



We used the stellar evolution code that is described in |Yoon, et| 
al. | (|2006) and references therein. It solves a set of stellar struc- 
ture equations using the Heney-type implicit method. The accel- 
eration term in the momentum equation was ignored to suppress 
pulsations that often occur because of partial ionization of hy- 
drogen when the surface temperature becomes lower than about 
5000 K (see, for e xample, |^ x>n & Cantiell o |2010 i. The opaci- 
ties are taken from Iglesias & Rogers |( |1996| l for T > 10 4 K, and 
from Alexander & Ferguson | ( |1994[ ) for lower temperatures. The 
nuclear network includes 63 nuclear reactions, with more than 
30 isotopes ( |Heger et al. |2000] l. The NeNa and MgAl chains are 
not considere d in this study. We adopted the rate of Caughlan 

I I I n , . 1 £. _ .... J - 



et al. 



I985| l for the l2 C(a, y) ls O reaction, multiplied by a fac- 



tor of 1.6. The considered reactions for carbon burning include 
2 12 C + 7 -> 20 Ne + a, and 2 12 C + y -> 23 Na + p, with sec- 
ondary reactions like 23 Na(p, Q-) 20 Ne and 23 Na(p, 7) 24 Mg. Neon 
burning is mainly followed by the reactions 20 Ne(7, or) 16 0, and 



20 Ne(o', 7) 24 Mg. For oxygen burning, only 28 Si and 32 S were 
considered as the main products. The nuclear network is not cou- 
pled with the transport equation. This may have consequences 
in the stellar models where protons are mixed into the convec- 
tive layers associated with helium shell burning, as discussed in 
Sect.g] 

The effect of rotation on the stellar structure is considered 
following Endal & Sofia~|( [r976| l. The transport of angular mo- 
mentum and chemical species through rotationally induced hy- 
drodynamic instabilities are considered as explained in Heger 
|et al. | ( [20 00). The effect of magnetic torques is also consid- 
ered in our models. In the radiative layers, dynamo action can be 
achieved by the interplay between amplification of toroidal fields 
(B$) by differential rotation, and continuous creation of radial 
fields (B x ) by the Tayler instability (the so-called Spruit-Tayler 
dynamo; Sprui t |2002| l. To consider this effect, the steady-state 
solution for and B T given by Spruit was used to calculate the 
effective viscosity for the radial transport of angular momentum 



(see Petrovic et al. 2005 for more details on the implementation 



of the Spruit-Tayler dynamo in our code). Even though the ef- 
ficiency of the Spruit-Tayler dynamo is currently much debated 
(e.g.,|Braithwaite |2006[|Zahn et al. |2007| |Gellert et al. |2008 



|Arlt & Rudige r |2011| l, it still remains one of the most promising 
mechanisms that can explain the observed spin rates of young 
neutron stars and isolated white dwarfs, which should reflect the 
history of angular momentum redistribution in the evolution of 
their progenitors ( |Heger et aT~||2005| |Suijs et al.~||2008[ ). Note 
that very weak magnetic fields suffice for the seed of the Spruit- 
Tayler dynamo ( |Spruit |fl999| |2002| l. Since there exist various 
plausible mechanisms for the generation of the primordial mag- 
netic fields in the early Universe (see Kandus et al. ||20PT for 
a recent review), our assumption of the presence of weak seed 
magnetic fields in the first stars can be justified. 

For convective overshooting and rotationally induced chem- 
ical mixing, we adopted the parameters given by |Brott et al. 
( |201 l| l, who have recently calibrated the mixing efficiency us- 
ing the data of the VLT- FLAMES Survey of massive stars. We 
used 0.1 and 0.028 for the and f c parameters, respectively 
(see |Heger et al. |[2000| for the definition of these parameters). 
We applied overshooting for the hydrogen-burning core only, 
for which a parameter of 0.335//p was adopted. We ignored 
overshooting for the stars beyond the core hydrogen exhaustion, 
where the efficiency of overshooting is expected to be small. 
This is because the convective cores during the advanced stages 
tend to grow in mass, which creates strong compositional strat- 
ification at the upper boundary of the convective core (jLanger 
[T99TT >. 

2.2. Mass loss 



|Krticka and Ku bat (2006) found that mass loss due to stellar 
winds from metal-free hot massive stars is practically negligi- 
ble. A very low mass loss rate of 10~ 14 M Q yr _1 is predicted only 
when the star reaches the Eddington limit. Therefore, we consid- 
ered stellar wind mass loss only if there is the surface enrichment 
in CNO elements, as the following. For hydrogen-rich stars, for 
which the surface hydrogen mass fraction (X S (H)) is greater than 
0.3, we used th e mass loss rates of Kudritzki e t al. I ( fT989| ) if 



r s > 10 K and Ni euwenhuijzen & de Jager ( 1990! other wise, 
with a metallicity dependence of Z u oy . On the other hand, Vink 
|& de Koter~| ( |2005[ ) found that the CNO lines play the dominant 
role for driving winds from helium-rich stars (i.e., X S (H) < 0.3) 
forZ < 10 3 Z , while the iron lines become important for higher 
metallicities. In our models, therefore, the mass loss enhance- 
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Fig. 1. Evolution of the models in the HR diagram for 10, 15, 20, 30, 60 and 100 M Q . Evolutionary tracks for different initial 
rotational velocities at a given initial mass are given in each panel as indicated by the labels. The end points of core hydrogen 
burning and core helium burning are marked by an square and a cross, respectively, on each evolutionary track. 



ment caused by the surface enrichment of CNO elements for 
helium-rich stars is considered only up to Z = 10 s , following 
the prescription by Nugis & Lamers (2000), which has a metal - 
licity dependence of Z 3 . 

Some cautionary remarks should be made here. The recent 
studies by |Krticka and Kubat] ( |2009| l and |Muijres et al. ipOTT) 
show that the effect of the surface enrichment of CNO elements 
on the mass loss enhancement from very metal-poor stars is 
only moderate. This means that the above prescription on the 
radiation-driven wind may overestimate the mass loss rates in 
our metal-free massive star models. In particular, CNO elements 
are highly ionized for very hot stars with surface temperatures 



higher than 50000 K, providing too few efficient lines for driv- 
ing winds. As we will see below, the helium-rich stars resulting 
from the chemically homogeneous evolution in our model grid 
have surface temperatures significantly higher than 50000 K. 

However, this uncertainty should not have any significant 
consequence for the conclusions of the present study, because 
the mass loss history in all our rotating model sequences is dom- 
inated by the mechanical winds that occur at the so-called a>T 
limit, which is the critical limit where the effective gravity at 
the stellar surface becomes zero because of the radiation pres- 
sure and the centrifugal force ( Langer 1 1997 ; Maeder & Meynet 
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Fig. 2. Same as in figure[l] but for 150, 200, 250, 300, 500, and 1000 



2000a). This effect is treated as in the following equation: 



/ 1 \ ' 43 

M(v) = max[l0~ 14 M G yr',M(v = 0)]l- — - 



where 



v GM 
Q := and v crit = J -^-(1 - O 

Vrrit V R 



(1) 



(2) 



Here, v denotes the rotational velocity at the equatorial surface, 
M(v = 0) the mass loss rate given by our adopted prescription for 
non-rotating stars that is described above, and T the Eddington 
factor. 

The enhancement of the mass loss rate when the star ap- 
proaches the (liF limit is reflected in the term (1/(1 - O)) 43 . We 



evaluate F using the radiative opacity for the outermost layers 
where the optical depth is lower than 100. Because a singularity 
can occur at CI = 1 according to Eq. [T] we limit the mass loss 
rate such that the mass loss time scale is longer than the thermal 
time scale of the star as in | Yoon et al. | ( |20T0) l: 



M = min 



M(v),0.3- 



M 



(3) 



where tkh is the Kelvin-Helmoltz time scale, and M(v) is given 
byEq.[T| 

The wind material that is mechanically ejected at the air 
limit may form a decretion disk. Recently, Krticka et al. ( 201 \) 
argued that viscous coupling between the star and the decretion 
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disk may significantly reduce the mass loss rate. This might po- 
tentially alter the final masses of the rotating models, but its de- 
tailed investigation is beyond the scope of the present study. 

2.3. Initial conditions and model grid 

The initial conditions and some important properties of the evo- 
lutionary models are summarized in Tables [T] and [2] We con- 
sider twelve different initial masses (M; n i t = 10, 15, 20, 30, 60, 
100, 150, 200, 250, 300, 500, and 1000 M ). Several different 
initial rotational velocities were considered at each mass, from 
zero to 80% of the Keplerian value (vk) at the equatorial sur- 
face. Solid-body rotation was assumed for the initial velocity 
profile. The initial mass and the initial rotational velocity are de- 
noted by the model sequence number. For example, the sequence 
"ml00vk04" indicates M init = 100 M Q and v illit /v K = 0.4. 

In metal-free massive stars, the CNO cycle cannot be ac- 
tivated initially. Because the energy generation due to the pp 
chain is too weak to support a massive star with M > 20 M Q 
for a significant fraction of the evolutionary time, the stellar core 
rapidly contracts until enough carbon (i.e., X(C) ~ 10 -10 ) is pro- 
duced by helium burning at T c ~ 10 8 K. Hydrogen burning with 
the CNO cycle only begins thereafter. This makes it difficult to 
make a relaxed model on the zero-age main sequence for Pop 
III stars. For this reason, we included a small amount of 3 He for 
the initial chemical composition, which may result from deu- 
terium burning, such that X(H) = 0.76, X( 4 He) = 0.23999 and 
X( 3 He) - 0.00001. In this way, all our initially relaxed models 
are supported by the energy generation due to 3 He burning, and 
the radii of our initial models are much larger than those obtained 
when the main sequence begins. Although this initial configura- 
tion looks artificial, it is not unrealistic given that the first stars 
should undergo rapid mass accretion phase that leads to signif- 
icant expansion of the envelope before they reach the zero-age 
main sequence (ZAMS) ( |Omukai & Palla |2003||Ohkubo et al~| 
|2009] l. 

In Table [T] we provide the information for the mass and ro- 
tational velocity when the stars arrive on the ZAMS, as well as 
their initial values. Given the efficient transport of angular mo- 
mentum by the Spruit-Tayler dynamo, near solid-body rotation 
is maintained during this initial contraction phase. It should be 
noted that Ekstrom et al. (2008 ) determined rotational velocities 
only when the stars are close to the ZAMS, and thus our ZAMS 
values should be compared with theirs, and not with our initial 
values. 



3. Basic results 

3.1. Evolution on the Hertzsprung-Russell diagram 

Figs. [T] and [2] show the evolutionary tracks of all of our mod- 
els in the Hertzsprung-Russell (HR) diagram. As mentioned 
above, our initially relaxed models are supported by 3 He burn- 
ing. Once the evolution starts, 3 He is very quickly destroyed and 
the stars undergo rapid contraction until they reach the ZAMS. 
The time scale of this initial contraction phase ranges from 10 4 
to 3 x 10 5 yr, depending on the initial mass. The models with 
Mnit = 10 and 15 M are powered by pp-chains for a significant 
fraction of the main sequence lifetime (13.3 Myr and 2 Myr, re- 
spectively, for non-rotating models). Thereafter, the CNO cycle 
becomes dominant, which is marked by the smooth turn-over to- 
ward lower temperature in the HR diagram. For M mlt > 20 M Q , 
the stars are already supported by the CNO cycle on the ZAMS, 



which agrees well with the previous work (Marigo et al. 2001 ; 
|Ekstrom et al. |2008| l. 

Compared to the models in the previous work, our non- 
rotating models become redder in the HR diagram during the 
post-main sequence phases. As an example, the surface temper- 
ature of the non-rotating 15 M Q models becomes log T e ff - 4.3 
and 4.5 in Marigo et al. and Ekstrom et al., respectively, while 
it becomes as low as log T e g = 3.9 in our case. In addition, our 
non-rotating models with M > 150 M Q become red-supergiant 
stars already during core hydrogen burning, while M > 250 M Q 
was needed in Marigo et al. for this to happen. This must result 
mainly from the fact that our adopted overshooting parameter 
is larger: ~ 0.3 Hp and 0.2 Hp are used in Marigo et al. and 
Ekstrom et al., respectively, while 0.335 Hp in the present study. 
The larger overshooting parameter also has an impact on the size 
of the stellar core and the evolutionary time. For example, the 
main sequence lifetime (th) and the helium core size (Mee-Core) 
of our non-rotating 200 M G model are 2.40 Myr and 1 14 M , 
respectively, (see Table [2]), compared to 2.27 Myr and 102 M 
in Ekstrom et al. 

Rotating models start at lower surface temperature and lu- 
minosity in the HR diagram for larger Vi n i t /vK, reflecting the ef- 
fect of the centrifugal force on the stellar structure. Rotationally 
induced mixing, on the other hand, tends to increase the size 
of the convective core and the surface helium abundance. The 
latter effect becomes dominant as the star evolves, leading to 
higher luminosity and higher surface temperature than those in 
the corresponding non-rotating models, for most of the time on 
the main sequence. The evolutionary time also becomes system- 
atically longer for higher initial rotational velocity because of the 
larger convective core. The central temperature is initially lower 
for higher initial rotational velocity because of the effect of the 
centrifugal force, but becomes higher already in the early phase 
of core hydrogen burning than those in a corresponding lower 
velocity model, as the convective core becomes larger in the ro- 
tating models due to rotationally induced chemical mixing. For 
instance, the center of the 100 M Q star enters the pair instability 
regime for Vj n i t /vK = 0.6, while it does not in the corresponding 
non-rotating sequence (see Fig. [3]). The final fate of such a mas- 
sive star can be significantly influenced in this way, as discussed 
in Sect. [6] 

One of the most dramatic effects of rotationally induced 
chemical mixing may be the so-called chemically homogeneous 
evolution (CHE; [Maeder | [1989), which can be observed with 
the blueward evolutionary tracks in Figures [T] and [2] We find 
that CHE can occur for a certain range of the initial mass 
(15 < Mi n it < 200 M Q ) if the initial rotational velocity is suffi- 
ciently high. For some cases, the star follows the CHE track for 
most of the main sequence phase but afterwards moves rapidly 
to the right in the HR diagram, in particular during the post- 
main sequence phases. Such transitionary evolution is found in 
Seqs. m200vk04, m250vk04 and m250vk06, for example. In the 
present study, we define these three different evolutionary pat- 
terns, i.e., normal evolution (NE), CHE, and transitionary evolu- 
tion (TE), according to the following criteria. 

- NE: The mass fraction of helium at the surface (T s ) remains 
smaller than 0.7 until the end of the main sequence. The star 
generally moves redwards in the HR diagram throughout the 
evolution. 

- CHE: T s becomes larger than 0.8 by the end of the main se- 
quence. The star generally moves bluewards throughout the 
evolution. 
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model sequences for Vinit/VR = (upper panel) and Vinit/vic = 0.6 
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- TE: F s becomes larger than 0.7 by the end of the main se- 
quence. The star follows the CHE track until the end of the 
main sequence, but moves redwards during the post-main se- 
quence phases. 

We discuss the conditions for CHE in more detail below 
(Sect. [5). 

3.2. Internal structure 



Fig. A.l show the internal structure for model sequences with 
Vinit/vK = 0.0, 0.2, 0.4 and 0.60 Note that the relative size of the 
convective core compared to the total mass becomes larger for 
higher initial mass, which is the well-known effect of higher ra- 



1 In some model sequences, the convective layers in the hydrogen 
envelope during the late evolutionary stages look very noisy (m20vk02, 
m60vk02, and ml5vk04). This is because the onset of convection de- 
pends sensitively on a minute change of the chemical gradient that is 
present in these layers. This phenomenon is frequently observed in mas- 
sive star models. 



diation pressure in more massive stars. We find that the convec- 
tion zone in the envelope extends up to the near-surface layers 
during the post-main sequence phases in all of our non-rotating 
models except for M ln \ t =15 and 20 M Q . This outer convec- 
tive zone is not developed for M < 60 M Q in Ekstrom et al. 



(2008 1 for their non-rotating case. This difference must originate 
from the larger overshooting parameter adopted in the present 
study than in Ekstrom et al., as explained above. The convec- 
tive carbon-burning core is not developed for M > 20 M in our 
models, which agrees well with the previous studies (Chieffi & 
Limongi |2004{ |Ekstrom et al. |2008[|Heger et al. |2000| i. 

Interestingly, in non-rotating models of M > 200 M G , the in- 
ner boundary of the convective envelope begins to penetrate into 
the layers that have previously undergone hydrogen shell burn- 
ing during core helium-burning phase. This eventually leads to a 
dredge-up of helium core material, and hot bottom burning oc- 
curs afterwards where primary nitrogen is abundantly produced. 
We discuss this phenomenon below in greater detail (Sect.|4jl. 

Rotating models have a larger helium core and a smaller con- 
vective envelope than those in the corresponding non-rotating 
models, which is an effect of rotationally induced chemical mix- 
ing. Dredge-up of helium core material by the downward pene- 
tration of the convective hydrogen envelope does not occur any 
more with rotation. The envelope remains radiative in the CHE 
models. 



3.3. Rotation and angular momentum 

Figure [4] shows the evolution of the surface rotational velocity 
in the rotating model sequences for 30, 100, 250 and 500 M . 
It rapidly increases during the initial contraction phase, and de- 
creases for a short period when the star begins to expand on the 
ZAMS. Note that in model sequences with Vj n ; t /v K > 0.4, the 
critical rotation is achieved during this initial contraction phase, 
leading to mass shedding even before the star arrives on the 
ZAMS. This is the reason why, in Table [T| the masses on the 
ZAMS for these rapidly rotating models are somewhat smaller 
than the initial masses. 

Soon after the ZAMS, the surface rotational velocity gradu- 
ally increases again, which is an effect of the internal transport 
of angular momentum. On the main sequence, the density of the 
hydrogen-buming core increases as the mean molecular weight 
increases as a result of nuclear burning. Although this tends to 
create differential rotation between the core and the envelope, 
magnetic torques due to the Spruit-Tayler dynamo keep the star 
rotating almost rigidly by transporting angular momentum from 
the core to the envelope. The outer layers are spun up accord- 
ingly, and most of our rotating models eventually reach the crit- 
ical rotation on the main sequence. The critical rotation at the 
surface is achieved earlier and maintained for a larger fraction 
of the evolutionary time, for a higher initial rotational velocity. 
If the star becomes a red supergiant afterwards, the expansion of 
the envelope leads to a rapid decrease of the surface rotational 
velocity (v/v crit < 0.1) . With CHE, on the other hand, it still 
remains close to the critical limit (v/v cl -i t > 0.8) during the post- 
main sequence phases. 

As prescribed in Eqs.[T[|2]and[3] our models lose mass when 
the surface velocity approaches the critical limit. As shown in 
Tables[TJand|2] more mass is lost for higher initial rotational ve- 
locity. In the case of M ml( = 250 M Q , for example, the total mass 
ejected (AM e j) is only about 1.0 M Q for v'i n ; t /vK = 0.2 (Seq. 
m250vk02), while it increases to about 60 M for v; n i t /v K = 0.6 
(Seq. m250vk06). Our rotating models generally lose more mass 
than those of Ekstr om et al. | ( |2008[ ) for a given initial rotational 
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velocity. For instance, our 60 M star in Seq. m60vk02 has a 
lower surface rotational velocity of 631 km s on the ZAMS 
(see Table [I} than that of in Ekstrom et al. (800 km s _1 ), but the 
total amount of the ejected mass is significantly larger (AM e j = 
3.58 M ) than that of Ekstrom et al. (AM ej = 2.32 M ). This 
is because our models reach the critical rotation earlier due to 
the Spruit-Tayler dynamo than those of Ekstrom et al., who did 
not include the effect of magnetic fields. Mari go et al~| (2003), 
who assumed rigid body rotation in their models, adopted a 
low initial rotational velocity: 500 km s -1 , which roughly cor- 
responds to vzams/vk ^~ 0-1 according to the figures 6 and 7 
in their paper. Their models therefore reach the critical rotation 
only briefly near the end of the main sequence, losing much less 
mass than our rotating models for a given initial mass. For exam- 
ple, the 250 M model in Marigo et al. gives AM ej = 3.35 M , 
while we have AM e j = 16.4 M in Seq. m250vk02, for which 
v'zams/v'k = 0.36 (TablefT). 

It is believed that the final fate of a massive star can be 



critically influenced by rotation of the core (e.g. Heger et al. 
2000). For example, the so-called collapsar (Woosley | |1993| l is 
predicted to occur if the final angular momentum in the core 
is higher than the critical limit for the formation of an accre- 
tion disk around a rotating black hole at a given mass (jKen-.iso; 
see also the figure caption of Fig. [5] |Bardeen et al. ||1972|l. 



m200vk06, m500vk04). In Seq. m30vk02, where the star gen- 
erally expands throughout the evolution, the core experiences 
strong braking due to the core-envelope coupling via the mag- 
netic torque. The final angular momentum in the core becomes 
well below y'Kerr.iso- In Seq. m30vk06, on the other hands, the 
star of the same initial mass undergoes CHE, to become a mas- 
sive helium star by the end of the main sequence and to avoid 
the core-envelope couping. The specific angular momentum in 
the core remains higher than jKerr.iso until the end of calcula- 
tion (core neon exhaustion), meaning that this model meets all 
the requirements for long GRB progenitors within the collapsar 
scenario (Woosley 1993): high angular momentum in the core, 
large core mass to form a black hole, and absence of an extended 
hydrogen envelope. This implies that the CHE channel for the 
formation of long GRB progenitors can still be important for the 
first stars (see Yoon & Langer 2005 1 Yoon, et al. 2006 1 Woosley 



& Heger 20061 for a more detailed discussion on GRB progen 



Figure|5]shows the evolution of specific angular momentum pro- 
file (j' r ) in four different model sequences (m30vk02, m30vk05, 



itors via CHE). Note that the condition for mag netar formation 
(j > 4 x 10 15 cm 2 s _1 ; e.g. jWheeler et al. |2 000) is also satisfied 
in this model. 

However, CHE does not always lead to retention of a large 
amount of angular momentum in the core such that j > y'Kerr.iso- 
In Seq. m200vk06, the specific angular momentum in the core 
decreases significantly below jKerr.iso before core helium exhaus- 
tion. Given the final mass of 146.43 M , the star of this model 
sequence is supposed to explode as a pair instability supernova, 
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Fig. 7. Ratio of the total number of ionizing photons from non- 
rotating models to that from rotating models with v^t/vx = 0.4. 
The connecting lines of filled circles, triangles and squares give 
the values for hydrogen and first and second helium ionization, 
respectively. 



but even if collapse into a black hole occurred, a GRB would 
not be produced in this case. We find that the CHE models have 
lower specific angular momentum in the core for larger initial 
mass, in general. This can be understood by the fact that the stel- 
lar structure of a rigidly rotating chemically homogeneous star 
is such that j r is smaller at a given mass coordinate, for a given 
v/ vk as shown in Fig. [6] In the CHE models, the surface ro- 
tational velocity remains close to the critical rotation for most 
of the lifetimes (Fig. HI, while near-rigid rotation is maintained 
until core helium exhaustion. Significant differential rotation is 
developed only thereafter. In addition, higher mass stars have 
lower v/vk at the critical rotation because of the increased role 
of radiation pressure. The general tendency for smaller j r with 
larger mass in the CHE models results from these factors. 

For Mjnit > 250 M Q , CHE is not realized anymore, and ex- 
pansion of the envelope to the red supergiant phase leads to 
strong core-braking as shown Fig. [5] As a result, the so-called 
super-collapsar (collapsar in stars of M > 250 M Q ) is not likely 
to occur. We discuss below the implications of our results for the 
final fate of the first stars in greater detail (Sect. [6]). 

3.4. Ionizing flux 

Detailed information on the hydrogen and helium-ionizing flux 
emitted from our models is provided in Table[3] Black-body radi- 
ation is assumed for the calculation of energy and number of ion- 
izing photons. Compared to Heger & Woosley (2010), our non- 
rotating models at a given mass produce more ionizing photons, 
which results from the larger overshooting parameter adopted in 
our study. Rotation tends to make stars brighter, bluer and live 
longer, and thus more ionizing photons are emitted than in the 
corresponding non-rotating case, in general. In Fig. [7] the to- 
tal number of ionizing photons from the non-rotating models is 
compared to that from the rotating models with Vi n ; t /vK = 0.4, 
for hydrogen (H), neutral helium (He) and singly ionized helium 
(He + ). The most notable change with rotation is found for ion- 
ization of He + . Here, the models with 20 - 150 M follow CHE, 
and the number of He + ionizing photons increases by factors of 



4 - 20. For H and He, enhancement of the number of ionizing 
photons is less prominent, but it is still significant for the CHE 
models, showing factor of 2 - 4 increase. 

4. Effects of chemical mixing on the 
nucleosynthesis 

The heavy elements synthesised in Pop III stars can be ejected 
into the surrounding medium by stellar winds and/or supernova 
explosions, which is supposed to critically influence the chem- 
ical evolution of the early Universe. Tables HJ and [5] present 
the yields for eight different isotopes in the wind material that 
has been ejected from the star throughout the evolution and in 
the last computed model, respectively, in each model sequence. 
Our non-rotating models do not experience mass loss at all, and 
therefore the wind yields of non-rotating models are not given 
in Table BJ The model sequences with M; n i t = 1000 M Q are 
not included in the tables either, because the calculations were 
terminated before the end of main sequence. The end points of 
the other models are mostly carbon burning or carbon exhaus- 
tion (see Table [2|, and none of our models have undergone the 
late evolutionary phases beyond oxygen burning. Therefore, in 
Table|5] neither compact remnant mass nor explosive nucleosyn- 
thesis is considered. This means that the amounts of 12 C, 16 
and 22 Ne in Table[5]should be taken as an upper limit for the su- 
pernova yields since they would decrease during the latest evo- 
lutionary stages and the supernova explosion that are not fol- 
lowed in this study, while the yields for the other isotopes that are 
mostly located in the outer layers of the star would not change 
much from the end points of our calculation. 

In Table [5] we also give the neutron excess at the center of 
the star (77,,) in the last model. Here, the neutron excess is defined 
as 

17 = Yj(Ni - Zi)Yi , (4) 

where N,, and F, denote neutron number, proton number, and 
abundance of isotope i, respectively. The last column of the table 
gives the maximum mass fraction of 14 N ( X max ( 14 N)) that has 
been achieved during the evolution. 
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In the following subsections, we discuss how convection and 
rotationally induced chemical mixing influence the nucleosyn- 
thesis in our models. 



4. 1 . Convective mixing 

During the post-main sequence phases, convection can induce 
mixing of helium -burning products (notably 12 C and 16 0) into 
hydrogen-shell-burning layers, boosting the CNO cycle. We find 
that there exist two different ways for this mixing. The first one is 
penetration of the helium convection zone (either convective he- 
lium core or convective helium-burning shell) into the hydrogen- 
burning shell, which has also been found by Ekstrom et al. 
(2008| and Heger & Woosley (2010 1. The second way is pene- 
tration of the convective hydrogen envelope into the helium core, 
which occurs for Mi n j t > 200 M Q in our models. To our knowl- 
edge, the latter has not been reported in the previous work. 

Penetration of the convective helium burning shell into the 
hydrogen-burning shell occurs for Seqs. m20vk02, m30vk00, 
mlOOvkOO and ml50vk02. In the Kippenhahn diagram of Seq. 
m20vk02 (Fig. | A. 1 [ >, for instance, the convection zone of the 
helium-burning shell continuously extends upwards and reaches 
the bottom of the hydrogen-burning shell when ff — t 10 yr. The 
consequent mixing of carbon into the hydrogen-burning shell in- 



duces boosting of the CNO cycle, and the hydrogen shell be- 
comes convective. This produces a large amount of primary ni- 
trogen. X max ( 14 N) in Seq. m20vk02 becomes as large as 10~ 3 
(see Table [5j, which is 1000 times larger than in the case with- 
out this mixing (i.e., X max ( 14 N)^ 10~ 6 in Seqs. m20vk00 and 
m20vk03). In terms of the total yield of nitrogen (Mu N ) it is 
about 10000 times higher than in Seqs. m20vk00 and m20vk03. 
However, M» N remains smaller than those in the corresponding 
CHE models that experience strong rotationally induced chem- 
ical mixing: e.g., M« N st 10~ 3 M Q in Seq. m20vk02, while 
Mw N a 10~ 2 M G including the wind yield (Tables [i] and El. 

In these sequences, a small amount of hydrogen is also mixed 
into the convective region of the helium-buming shell. Given the 
large mass fraction of carbon and the fairly high temperature 
(~ 2 x 10 8 K) in this region, the nuclear time scale for hydro- 
gen burning becomes shorter than the convective turn-over time 
scale. This rapid nuclear burning in convective layers cannot be 
properly described with our code because the nuclear network 
and the transport equation are calculated in succession. This 
leads to an overestimate of the energy generation rate here since 
the hydrogen mass fraction at the bottom of the convection zone 
above the helium burning shell can become artificially higher 
than it should be. However, we find that the overall amount of 
nitrogen and energy produced in this region is small compared to 
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that of the hydrogen-burning shell. In addition, our calculations 
are terminated soon after this mixing occurs, and this uncertainty 
should not affect the main conclusions of the present study. 

Penetration of the convective helium-core into the hydrogen- 
burning shell can be most prominently observed in Seqs. 
mlOOvkOO and ml50vk02. Here, the convective helium core 
grows and reaches the bottom of the hydrogen-burning shell, as 
shown in Fig. |A.l| In this case, a significant fraction of the newly 
produced primary nitrogen in the hydrogen shell source is mixed 
into the convective helium core, enhancing the neutron excess in 
the core, which is discussed below in Sect. |4.3| 

On the other hand, our non-rotating models with M; n ; t > 
200 M experience dredge-up of the helium core material into 
the hydrogen-burning shell by convection in the hydrogen enve- 
lope (see Fig. |A. 1[ >. As an example, the detailed structure around 
the boundary between the helium core and the hydrogen enve- 
lope is shown in Fig. [8] for Seq. m500vk00. The hydrogen en- 
velope becomes convective already on the main sequence and 
the convection zone progressively extends downwards. During 
core helium burning, it eventually penetrates the layers that pre- 
viously belonged to the convective helium core. A large amount 
of carbon and oxygen, which has been produced by core helium 
burning, is thus mixed into the hydrogen envelope, boosting the 
CNO cycle at the bottom of the convection zone. Fig. [9] shows 
that the mass fraction of the CNO elements in the 500 star model 
becomes as large as 0.03, after this dredge-up history. 

We find that this convective dredge-up of the helium core 
material becomes more efficient for more massive stars, in gen- 
eral. The total yields of 13 C and 14 N increase from 1.5 x 10~ 3 M Q 
and 3.71 x 10~ 4 M in 200 M Q model to 0. 0.34 M Q and 
0.12 M in 500 M Q model, respectively. Note that in our non- 
rotating models, the outermost layers of about 0.1 M Q remain 
radiative even when the most part of the hydrogen envelope be- 
comes convective. Since chemical mixing cannot occur beyond 
the convection zone without rotation, the initial values are main- 
tained for the chemical composition at the surface. This is the 
reason for the sudden change in the mass fractions of the iso- 
topes near the surface in Fig. [9] No mass loss occurs in our 
non-rotating model sequences, given our wind prescription that 
allows stellar winds only with surface-enrichment of heavy ele- 
ments (Sect.|2~2]i. In reality, massive stars may become unstable 



to the pulsational instability, which may cause mass loss ( Baraffe 
et al. 2005 ). Pulsation caused by partial ionization of hydrogen 



may be particularly important during the red supergiant phase, 
when the above-discussed convective dredge-up can occur ef- 
ficiently. Therefore, we cannot exclude the possibility of sig- 
nificant mass loss from non-rotating massive Pop III stars, and 
the consequent enrichment of the CNO elements into the cir- 
cumstellar material. Because pulsation-driven winds from RSG 
stars must be slow, it should have interesting consequences in 
the chemical evolution history as in the case for massive AGB 
stars (e.g., Ventura et al. |200T) . Even formation of dusts around 
such Pop III RSG stars might be possible, given the high abun- 
dance of CNO elements in the envelope (Fig. |9j. This question 
deserves more studies. 



4.2. Rotationally induced mixing 

The significance of rotationally induced chemical mixing on nu- 
cleosynthesis varies according to the initial rotational velocity 
for a given initial mass. As shown in Fig.[l0](see also Tables [4] 
and|5]l, in most of the slowly rotating models that follow NE, the 
total yield of primary nitrogen is not higher by more than a factor 
of 10 than in the corresponding non-rotating case, if the star does 
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not experience the encroachment of the helium convection zone 
into the hydrogen shell source. For instance, Seq. mlOOvkOO 
gives Mm n = 2.87 x 10~ 6 M Q , while Seq. ml00vk02 gives 
Mi4 N = 3.29 x 10~ 6 M . For models with M m = 30, 200, 
300 and 500 M , M» N is even larger in non-rotating sequences 
than in the corresponding sequences with Vi n j t /vK = 0.2. This 
is because these non-rotating massive stars undergo the con- 
vective mixing discussed above: penetration of the convection 
zone of the helium-burning shell into the hydrogen-burning shell 
for Mi n i t = 30 M , and the convective dredge-up of the helium 
core material for M init = 200, 300 and 500 M Q . On the other 
hand, Seqs. m20vk02 and m60vk02 give exceptionally large 
Mw N (8.14 x 10~ 4 and 9.21 x 10~ 5 M Q , respectively), compared 
to the other cases with the same initial velocity. This is not be- 
cause of the rotational mixing, but because of the convective 
mixing between the convective helium shell and the hydrogen 
shell source. 

The CHE models produce more primary nitrogen by several 
orders of magnitude than in the corresponding non-rotating or 
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slowly rotating case, in most cases. For instance, Seq. m30vk04 
gives Mh n = 0.024 M Q , compared to 1.19 x 10~ 7 M in 
Seq. m30vk00. Although they evolve almost chemically homo- 
geneously throughout the main sequence, some amount of hy- 
drogen is still left in the outermost layers, which undergoes hy- 
drogen shell burning during the core helium-burning phase (See 
Fig. |A.l| and[TT). Chemical mixing induced by rotation between 
the helium core and the hydrogen shell then leads to abundant 
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production of primary l4 N and 13 C, as well as 22 Ne (see Fig 
for such an example). 

The total yields of primary nitrogen in the rotating models 
by Ekstrom et al. ( 2008b are compared with ours in Fig 
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In their study, where magnetic torques are not considered for 
the transport of angular momentum, none of the models fol- 
lows CHE, and the production of primary nitrogen is enhanced 
mainly by chemical mixing due to the shear instability during 
the core helium-burning phase, where a strong degree of dif- 
ferential rotation is developed across the boundary between the 
helium core and the hydrogen envelope. In our models with the 
Spruit-Tayler dynamo, the role of the shear instability is neg- 
ligible because near rigid-body rotation is maintained until the 
end of core helium burning. In our NE models, chemical mixing 
by the Eddington-Sweet (ES) circulations becomes rather inef- 
ficient during the post-main sequence phase, because the strong 
chemical stratification between the helium-burning core and the 
hydrogen envelope significantly slows down the ES circulations. 
This explains the weak production of primary nitrogen in our ro- 
tating NE models compared with the rotating models by Ekstrom 
et al. In the CHE models, however, the inhibition of the ES cir- 
culations by the chemical stratification is much less significant 
and the total yields of primary nitrogen are comparable to those 
of Ekstrom et al. The yields of other isotopes like 13 C, ls O and 
22 Ne also increase by several orders of magnitude with CHE. 

As shown in the Kippenhahn diagrams in Fig. A.l the CHE 
models with Mj„j t < 150 M Q lose the hydrogen-burning shell 
layers before the end of the core helium-burning phase. Note that 
the mass loss here is dominated by the mechanical winds, which 
must be slow. Slow ejecta of hydrogen-burning products via ei- 
ther winds from massive AGB stars or mechanical winds from 
massive single or binary stars are often invoked to explain some 
chemical anomalies observed in globular clusters, characterized 
by enrichment of He, N, Na and Al, and depletion of C and O 
(| Ventura et al. |2001[|Denissenkov & Herwig |2003||Prantzos & 
ICharbonnel |2006||Decressin et al. |2007Hde Mink et al. |2009> . 
As shown in Table |4] the wind materials of our CHE models are 
indeed marked by enrichment of He and N, but the amounts of 
C and O may be too large to explain the anomalies. We will ad- 
dress the implications of these models for the chemical evolution 
in the early Universe in a future paper. 
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Fig. 10. Total yields of nitrogen as a function of the initial mass 
for Vi n it/vK = 0.0 (star symbol), 0.2 (open circle) and 0.4 (filled 
circle). For Vi n ; t /vK = 0.4, the CHE models are marked with a 
yellow color, while the NE/TE models with a purple color. The 
connecting line with open triangles gives the nitrogen yields of 
the rotating models by Ekstrom et al. ( 2008 1. 
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Fig. 11. Mass fraction of different isotopes in Seq. m30vk04, 
when primary nitrogen is abundantly produced during core he- 
lium burning. 



4.3. Neutron excess 

The neutron excess (Eq. Q in the stellar core is one of the 
most important factors that determine the nucleosynthesis of a 
supernova explosion. This is particularly the case for the pair- 
instability supernova (PISN). Detailed calculations by |Heger~& 
Woosley (2002) showed remarkable deficiency of odd-charge 



elements in the yields of PISNe from Pop III stars, which results 
from very low neutron excess. Despite the theoretical expecta- 
tion of frequent PISN explosions in the early Universe, such a 
peculiar abundance pattern has not been found in the observed 
very metal-poor stars, which are believed to conserve the signa- 
tures of the nucleosynthesis in the first generations of stars. It 
should be noted, however, that |Heger & Woosley~| (|2002) only 



calculated pure helium star models, ignoring the complicated 
history of chemical mixing that is discussed above. It is therefore 
important to investigate whether or not the convective and/or ro- 
tational mixing can significantly increase the neutron excess in 
Pop III PISN progenitors. 

In Table [5] we show the value of the central neutron ex- 
cess (77c) at the end of calculation for each evolutionary se- 
quence. In our calculations, the major source for the neutron 
excess is the I4 N(ff, y) 18 F(e + , v) I8 reaction, which can occur 
in helium-buming layers. This means that rj c in our models is 
determined by the amount of primary nitrogen mixed into the 
helium-burning core, and does not include any excess beyond 
helium burning. On the other hand, |Heger & Woosley"] ((2002 1 



include many more reactions that can further increase the neu 
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tron excess during the advanced stages beyond helium burning 
(e.g., the 20 Ne(p, y) 21 Na(e + , v) 21 Ne reaction during carbon burn- 
ing and many other weak processes). The implication is that the 
problem with the deficiency of odd-charged elements would not 
be solved by chemical mixing of primary nitrogen, unless the 
values of rj c in our models were significantly higher than those 
of Heger & Woosley. 

As discussed below (Sect. |6j, PISNe are expected for 120 < 
M init /M < 240 with NE, and 84 < M init /M < 200 with CHE. 
For this mass range, the neutron excess in our models is limited 
to ~ 10~ 5 , which is much lower tha n those of the Pop III star 
models by |Heger & Woosley | ( |2002"] ) (t] c ~ 10" 4 ). Therefore, we 
conclude that PISNe from our models would not result in a sig- 
nificantly different nucleosynthesis compared to that of Woosley 
& Heger's models. It should be investigated in future if this con- 
clusion would still hold with non-magnetic models, where strong 
mixing due to the shear instability would occur during the post- 
main sequence phases (Ekstro m et al. |20 08). 

5. Conditions for chemically homogeneous 
evolution 

The CHE has important consequences not only for the nucle- 
osynsthesis, but also for the final fate of massive Pop III stars 
(Sect. [6 1. We give the boundary for CHE in Fig. [T2]in the plane 
spanned by the initial mass and the surface rotational velocity 
|^|(see Sect. [6] for more detailed explanation on the figure). Our 
results indicate that CHE would not occur for Mzams ~ 1 3 M Q 
and Mzams Si 190 M . The limiting rotational velocity for CHE 
continuously decreases from Mzams = 13 M to about 60 M Q , 
and rapidly increases from about 150 M until it reaches the for- 
bidden region at Mzams — 190 M . The mass dependency of the 
limiting vzams for CHE can be understood as follows. 

In our models, chemical mixing is dominated by Eddington- 
Sweet (ES) circulations. According to Kippenhahn ( 1974), who 
we follow for the prescription of ES circulations in our models, 
the circulation velocity in a radiative region is given as 



(x> r 



V e ,K, 



where 

V'e,K : = 



L r (2(e B + e v )r 2 2r 2 



<5(V ad - V) GM r 



U 



M r 4npr 



Here a> r is the angular velocity at a radius r, w K ,r the local 
Keplerian value, e n the nuclear energy generation rate, e v the 
energy loss rate due to neutrino emission, and L r the local lu- 
minosity. The other symbols have their usual meaning. Note that 
v e becomes v 6j k when a> r = co^ r . 

The chemical stratification tends to slow down the circula- 
tion velocity, and a correction is made for the actual circulation 
velocity in our models ( Heger et al. |2000 1: 



v ES = max(|v e | - 1^1,0), 



where 



Hp 



(6) 



(7) 



TKH,r <J(V - V ad ) 

Here TKH.r is the local Kelvin-Helmoltz time-scale at a radius r. 



2 Note that in this figure, we use the values on the ZAMS for ro- 
tational velocity and mass (i.e., Vzams an d Mzams) given in Table [T| 
instead of the initial values. 



For the discussion below, we define three different time 
scales: 



tes.k := I —dM r IM md , 

J Mai Ve > K 

- f r 
T ES,0 •- I ~ 



■c/M r /M rad 



and 



J ,11 



-dM r /M md . 



(8) 



(9) 



(10) 



Here M rad denotes the total mass of the radiative layers in the 
envelope. Note that tes.k does not depend on the rotational ve- 
locity, and thus can be used as a "susceptibility" of a fluid to ES 
circulations, for a given thermodynamic condition. On the other 
hand, tes.o represents the ES circulation time scale for chemi- 
cally homogeneous layers since the effect of a chemical stratifi- 
cation is not included in v e . Given that tes.o ~ (wK, r /w r ) 2 TEs,K 
(see Eq. BJ, the ratio tes.o/tes.k gives information on how the 
ES circulation is influenced by rotation in stars. For example, 
loss of angular momentum and/or expansion of a star would re- 
sult in an increase of t E s,o/ t es,k- Finally, t E s is the real ES cir- 
culation time that considers the stabilizing effect of a chemical 
stratification. The ratio tes/tes,o can be used as a measure of 
chemical inhomogeneity inside a star: Tes/tes.o = 1 for = 
and tes / tes,o > 1 otherwise. 

It should be noted that the actual time scale of ES circula- 
tions in stellar evolution models is much larger than the above 
defined t E s, as a result of the calibration of the mixing efficiency 
with observational data (see Sect. [2}. Therefore, our discussion 
below should be considered as only indicative. 

We present the ratio of tes.k to the main sequence lifetime 
(? M s) as a function of initial mass in Fig. 13 to investigate how 



thermodynamic conditions vary according to the initial mass in 
favor/disfavor of ES circulations. Here, tes.k is estimated on the 
ZAMS for each mass. The figure shows that t E s,kAms is system- 
atically smaller for a higher mass star. This is related to higher 
radiation pressure and lower density in a higher mass star. This 
means that, for a given initial rotation velocity, CHE should be 
favored in a higher mass star. This may explain why the limiting 
rotation velocity for CHE continuously decreases with increas- 



ing mass from 13 M Q to about 60 M Q in Fig. 12 



However, the influence of mass loss becomes important for 

more massive stars 



(5) very massive stars: As discussed in Sect. 3.3 



can reach the critical rotation in an earlier evolutionary phase, 
losing more mass and angular momentum. This effect is re- 
flected in the evolution of Tes,o/i"es,k that is shown in the upper 
panel of Fig. [l4]for some model sequences with Vinit/VR = 0.4. 
Tes,o/ t es,k increases very rapidly in 250 M and 500 M Q stars, 
compared to those in lower mass stars. For example, in 500 M 
star, tes,o/t"es.k increases by a factor of two from its initial value, 
when the central helium mass fraction (Y c ) becomes about 0.4. 
Since the change of the radius during this period is only moder- 
ate (from 12 R Q to 15 R Q ), the rapid increase in tes,o/tes,k results 
mainly from the loss of angular momentum due to centrifugally- 
driven mass loss. 

The evolution of tes/tes,o is shown in the bottom panel of 
Fig. [14] For 30 M G and 100 M G models, t E s/t E s,o - 1.0 is 
maintained until the end of main sequence, as a result of CHE. 
For 250 M Q and 500 M , tes/tes,o starts deviating from 1.0 
at Y c =! 0.6 and Y c 0.4, respectively. This indicates that al- 
though these stars initially follow CHE, chemical mixing due to 
the ES circulation becomes inefficient later as they progressively 
lose large amounts of angular momentum. This interpretation 
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Fig. 12. Phase diagram for the final fates of massive Pop III stars, in the plane of the mass and the fraction of the Keplerian value 
of the equatorial rotational velocity on the zero-age main sequence. The dotted line denotes the limit above which the rotational 
velocity at the surface would exceed the critical rotation velocity (v cr it in Eq.|2j. The thick solid line marks the boundary between the 
regimes for the chemically homogeneous evolution (CHE), and the non-CHE evolution (normal evolution or transitionary evolution; 
see Sect. [3]). The regions for different final fates are divided by thin solid lines. Each acronym has the following meaning. SN IIP: 
Type IIP supernova, NS: neutron star, BH: black hole, SN II: type II supernova, PISN: pair-instability supernova, Puis. PISN: 
pulsational pair-instability supernova, GRB: gamma-ray bursts, HN: hypernova, SNIbc: type lb or Ic supernova. 



is confirmed by Fig. 15 which shows the evolution of surface 



helium mass fraction (Y s ) as a function of central helium mass 
fraction (T c ), for model sequences of m30vk04, m250vk04 and 
m500vk04. Note that T s is larger in 250 M G and 500 M Q stars un- 
til Y c =i 0.6 and 0.4, respectively, than in 30 M Q . This implies that 
chemical mixing is initially more efficient in those very massive 
stars than in 30 M Q star, although they eventually deviate from 
CHE because of angular momentum loss. 

This means that without angular momentum loss, the CHE 
regime in Fig. 12 would extend to larger masses. Specifically, 
(v k /vzams) 2 Tes,k, the result shown in Fig 
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given that r E s 

implies the limiting vzams/^k for CHE would decrease from 
0.74 for 20 M to about 0.4 for 500 M in Fig ~ 
no loss of angular momentum. 

From the above analysis, we draw the following conclusions. 



12 if there were 



- Thermodynamic conditions are generally more favorable for 
chemical mixing by the ES circulation in more massive stars. 
Therefore, the limiting rotational velocity for CHE can de- 
crease with increasing mass, as long as the ES circulation is 
not significantly slowed down through angular momentum 
loss. 

- Critical rotation can be more easily achieved in more massive 
stars, since their surface luminosity is close to the Eddington 
limit. The consequent rapid loss of mass and angular mo- 
mentum eventually tends to inhibit CHE in very massive 
stars (>~ 200 M ), even if rotationally induced chemical 
mixing is initially very efficient. 
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Fig. 13. t E s,k in Eq. 8] divided by the main sequence lifetime 
(?ms) and multiplied by 1000, as a function of the initial mass. 
The ZAMS models of 20, 30, 60, 100, 250 and 500 M were 
used for the evaluation of r E s K . 



6. Implications for the final fate 

In Fig. [12] we present a phase diagram of the expected final fates 
of rotating massive Pop III stars in the plane of the mass and the 
rotational velocity in units of the Keplerian value. Note that here 
we used the ZAMS values instead of the initial values for mass 
and rotational velocity (see Table T|and discussion in Sect. 2.3 1. 
To determine the boundaries for different regimes, we use inter- 
polated values between the grid points of our calculated models, 
when it was needed. 

The forbidden region in the figure is the regime where the 
surface rotational velocity exceeds the critical limit v cr j t (Eq.[2]l. 
Any star in this region would be unstable, experiencing mass 
eruption, and forced to come back to the regime of sub-critical 
rotation. The CHE regime (vche ^ ^zams ^ v crit) is marked by 
a thick solid line. Both NE and TE regimes are below this limit. 
Because our model grid is not fine enough to clearly identify the 
TE regime, we do not distinguish the two in the non-CHE regime 
in the figure. In any case, the parameter space for TE in this plane 
should be very small, located in the vicinity of the border of the 
CHE regime. 



6. 1 . Core collapse supernovae 

In the NE regime, the lower mass limit for black hole (BH) 
formation is simply assumed to be 25 M Q . We did not follow 
the evolution up to iron core formation, and we cannot discuss 
how rotation would change the pre-supernova structure, which 
may affect this limit. However, in all models in the NE regime, 
the amounts of angular momentum in the core fall below the 
limit for magnetar (i.e., j > 4 X 10 15 cm 2 s in the innermost 
1.4 M Q ; e.g., [Wheeler et al. |20Q0| > or collapsar formation (i.e., 
j > 4 x 10 16 cm z s" 1 in the layers that would be accreted to the 
black hole; Woosle y"|1993 i, by the end of core helium burning. 
This means that a jet-driven supernova or hypernova-like event 
powered by rotation may not occur in stars having hydrogen-rich 
envelopes, which have been discussed in Heger et al. ( 2003) 1 in 
detail. Type IIP supernovae via usual core-collapse to a neutron 
star are expected for 8 M Q < M Z ams ~ 25 M Q , and weak fall- 
back SNe or direct collapse to BH would occur for the other 
mass range outside the pair-instability regimes. 
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Fig. 14. Upper panel: The evolution of tes,o/t"es,k (see Eqs. [s] 
and[9]l in the model sequences m30vk03, ml00vk04, m250vk04 
and m500vk04. Only main sequence models are considered here. 
Bottom panel: The evolution of tes/tes,o ( see Eqs.[9]and 10 1 in 
the corresponding model sequences. 



The light curve, luminosity, and chemical yields of a SN 
IIP are critically affected by the radius of the progenitor star at 
the pre-supernova stage. Our 10 M models expand to the red- 
supergiant phase, having the final radii of Rf 450/?©. 15 and 20 
M models remain relatively blue up to core carbon exhaustion, 
having /? f - 130...150/? o and/? f =s 1 00 R Q , respectively. We find 
these values do not change much with rotation in the NE regime, 
but they may still be affected by uncertain factors such as core 
overshooting and production of primary nitrogen, as discussed 
in Sects. [3]andg] 

6.2. Pair-instability supernovae 

In the present study, we did not follow the explosive nuclear 
burning caused by pair instability and thus our models do not 
provide an exact boundary for the PISN regime. For the phase 
diagram of Fig. [12] we simply adopted the mass limits given by 
Heger & Woosley (2002): pulsational pair-instability SNe (puls- 
PISN) for helium core masses of 40 < M He coK /M Q < 64, and 
pair-instability SNe (PISN) for 64 < M H ^ COK /M Q < 133. Note 
that these limits are based on non-rotating stellar models, and 
that the effect of the centrifugal force can lead to an upward-shift 
of each boundary with rotation ( jGlatzel et al. |1985 1. Because in 
our potential PISN progenitor models the rotational velocity in 
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Fig. 15. Mass fraction of helium at the surface (Y s ) as a func- 
tion of the mass fraction of helium at the center (T c ) during the 
main sequence for model sequences m30vk04, m250vk04 and 
m500vk04. The thin dotted line corresponds to the case of per- 
fect chemical homogeneity. 



the innermost layers does not exceed a few percent of the local 
Keplerian values even in the case of CHE, the role of rotation 
may be minor. But more detailed studies should be carried out 
for clarifying the uncertainty^] 

In the NE regime, the mass limits for puls-PISN and PISN 
decrease with rotational velocity, because the helium core mass 
increases with rotation. In the CHE regime, they do not change 
much with rotation, since in this case, the whole star is trans- 
formed to a helium star by complete mixing by the end of the 
main sequence. The progenitors of PISNe in the NE regime have 
a very extended hydrogen envelope with Rf 3000... 5000 R Q , 
depending on the mass. As shown in the recent calculations by 
Kase n et al. | ( |201 1[ ), these PSNE would be marked by a long 
plateau phase of the light curve (100 - 400 days), and strong hy- 
drogen lines (absorption and/or emission) in the spectra. In the 
CHE regime, the progenitors of puls-PISNe and PISNe do not 
have any hydrogen. The total amount helium at the pre-SN stage 
is small, ranging from 0.3 to 0.7 M (Table [5}. A PISN from a 
CHE model is thus likely to appear as a SN 2007bi-like event, 
which is a very bright SN Ic ( |Gal-Yam et al. |2009| , if the final 
mass of the progenitor is high enough (> 90 M Q ) to produce a 
large amount of nickel (Heger & Woosley 2002; Kasen et al. | 
[20TT) . 



6.3. Gamma-ray bursts 

If specific angular momentum in the CO core is higher than 
7k,isoi an d if the helium core mass is not within the range of 
PISNe, a collapsar may occur. In the NE regime, no collapsar 
formation is expected because angular momentum is too low in 



Note added in proof: a paper by Chatzopoulos & Wheeler 1 2012 
has recently been posted on astro-ph. Their result on the mass range for 
pair instability supernovae from rotating Pop III stars agrees well with 
ours. 



the core, while the conditions for a collapsar are fulfilled for 
stars of 13 < M ZAMS /M Q < 84, in the CHE regime. The in- 
nermost cores (~ 1.4 M Q ) in the stars in the same regime also 
have high angular momentum to produce a magnetar (j >~ 
5.4 x 10 15 cm s ). Therefore, stars in this regime may be con- 
sidered as GRB progenitors within both the collapsar scenario 
( [Woosley |1993||Macfadyen & Woosley |1999|> and the magne- 
tar sc enario (e.g., IDuncan & Thomps on] 1992| Wheeler et al. 
2000[|Bucciantini et al. |2009| >r[ 

Interestingly, both puls-PISN and GRB are expected for 56 < 
A^zamsM^o $ 84 in the CHE regime. In this case, pulsational 
pair-instability that occurs after carbon burning in the core may 
eject several solar masses, probably episodically (Woosley et al. 
2007 ). The chemical composition of the puls-PISN ejecta would 
be characterized by large amounts of carbon and oxygen, and 
some helium (~ 0.3 M ). Only about several years thereafter, 
a GRB would be produced when the core collapses. Therefore, 
the GRB jet should interact with the very dense ejecta, which 
might result in an extremely luminous afterglow if the GRB jet 
were energetic enough. Heger et al. ([2003 ) also suggested the 



possibility of occurrence of both puls-PISN and GRB from the 
same progenitor for very low metallicity. However, in their sce- 
nario, the puls-PISN is supposed to eject the hydrogen-rich en- 
velope from the progenitor, while in our scenario based on the 
CHE models, no hydrogen would be contained in the ejecta of 
the puls-PISN. The accompanied supernova (if it were not out- 
shined by the GRB afterglow) would not appear as a SN 2006gy- 
like event (luminous type Iln Smith et al. |2007|> b ut rather look 
like a Quimby-type supernova (Quimby et al. 201 1 1 or a SN Ibn 
if helium were excited as a result of the interaction (Foley et al. 
|2"007l|Pastorello et al. 12007] ). 

The GRB progenitors with 12 < Mzams/M b < 56 may not 
suffer the pulsational pair-instability. In these models, however, 
the critical rotation is maintained at the surface for a significant 
fraction of the evolutionary time, which leads to significant mass 
loss (Fig. In particular, the mass loss rate becomes very high 
during the last 1 4 — 1 5 yrs, because of the rapid decrease of the 
radius of the star (Fig. [TJ. The mass loss would predominantly 
occur via a decretion disk in this case, and radiation driven winds 
along the poles must be very weak or negligible (see Sect. |2j. 
This may lead to a very complex, highly aspherical structure of 
the circumstellar medium, which might be imprinted on the af- 
terglow (see van Ma rie et al. |2008| for a related study). 

Recently several authors investigated possible observational 
signatures from collapsars in very massive stars with M > 
250 Mq (|Komissaroy & Barkov |2010]IMeszaros & Rees |2010l 
|Suwa & Ioka |2011[|Toma et al. |2011| i. Our result implies that 
a GRB-like event from such a super-collapsar may not occur. 
First of all, CHE is inhibited for very massive stars (Sect. [3J, 
and consequently they are expected to have a very extended hy- 
drogen envelope, from which a jet from the central engine may 
not break out easily. Even if CHE were realized, the final angu- 
lar momentum in such a massive star would not be high enough 
to produce a GRB (see discussion in Sect. 3.3 1. Suwa & Ioka 



(201 1 1 argued that for very massive stars, accretion of the long- 
lived hydrogen envelope would result in jet-breakout even if they 
had a supergiant envelope. However, our 250, 300 and 500 M 
models show that the specific angular momentum would become 



4 Explosions as an energetic supernova without making a GRB might 
be another possibility for a star in this regime, as implied by the numer- 
ical studies of B urrows et al. |p0 07 1 and Dessart et al. ( 2008 1 on mag- 
netically driven explosions caused by the magneto-rotational instability 
during core-collapse 
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well below v'k.Iso for the entire mass range including the hydro- 
gen envelope, at the pre-collapse stage (Fig. [5}. Direct collapse 
of the whole star into a BH without any remarkable visible event 
is the most likely outcome in this case. Although this conclusion 
is based on our assumption of the presence of weak-seed mag- 
netic fields in Pop III stars, our preliminary results show that 
necessary conditions for a super-collapsar still cannot be easily 
fulfilled even without the assumption of magnetic fields (Yoon 
et al., in preparation). 

According to theoretical studies on the formation of the first 
stars in the early Universe, Pop III stars may have two differ- 
ent classes in terms of initial mass function ^Johnson & Br omm 
[1 |2006l |Yoshida et al. ||2007l |Bromm et al. ||2009| l: the first- 
generation stars formed via collapse of neutral primordial gas in 
isolated dark-matter minihalos (Pop III.l), and the stars formed 
in an environment affected by photo-ionization from PopIII.l 
stars (Pop III.2). Formation of Pop III.l and Pop III. 2 stars is 
governed by different cooling agents (H2 and HD molecules, re- 
spectively) and Pop III. 2 stars are supposed to have much lower 
initial masses (~ 40 M Q ) than Pop III.l stars (~ 500 M Q ). In 
this framework, our result implies that the initial mass function 
of Pop III. 2 stars would be more favorable for GRB production 
than that of Pop III.l stars. 

7. Conclusions 



We presented a new grid of massive Pop III star models, which 
include the effects of rotation on the stellar structure and mass 
loss, rotationally induced chemical mixing, and the transport of 
angular momentum. The most up-to-date calibration for the mix- 
ing efficiency parameters by Brott et al. ( 201 1\ was adopted, 
and the presence of initial magnetic seed fields in Pop III stars 
was assumed such that the transport of angular momentum is 
dominated by magnetic torques due to the Spruit-Tayler dynamo 
(Spruit 2002). We summarize the main conclusions from the 
analyses of the model grid as follows. 



1 . Because of the larger overshooting parameter adopted here 
than in the previous work, our non-rotating models become 
redder on the HR diagram, and stars with M m \ t = 10 M 
and Mi n i t > 30 M become red-supergiants, even when they 
do not experience a boosting of the CNO cycle caused by 
mixing in the hydrogen shell source (Sect.[3jFigs.[T]and|2|. 

2. We find, for the first time, that dredge-up of helium core 
materials into the convective hydrogen envelope occurs 
as the bottom of the envelope moves down during core 
helium burning in non-rotating stars with Mj n i t > 200 M Q 
(Sect. |4j Fig. [BJ. This leads to mixing of large amounts of 
carbon and oxygen into the convective hydrogen envelope, 
and abundant production of primary nitrogen (Fig. |9]l . 
The stars are already in the red-supergiant phase when 
the dredge-up occurs. The hydrogen envelope at this stage 
must be unstable to pulsational instability driven by partial 
ionization of hydrogen. If pulsations caused mass loss at 
this stage, it would have interesting consequences such as 
enrichment of CNO elements and other hydrogen-burning 
products into the circumstellar region, and formation of 
dusts. 

3. The so-called chemically homogeneous evolution 
(CHE) is expected for a certain range of masses 
(13 < Mzams/M q < 190) and rotational velocities 
(Sects. [3] and |5j Fig. [12}. On one hand, thermodynamic 



conditions in more massive stars are more favorable for 
chemical mixing due to Eddington-Sweet circulations. On 
the other hand, more massive stars can reach the critical ro- 
tation at an earlier stage of the evolution to induce rapid loss 
of angular momentum, and Eddington Sweet circulations 
can be significantly slowed down. As a result, the limiting 
rotational velocity for CHE decreases from Mzams = 13 M G 
to about 60 M Q , but increases from ~ 150 M to ~190 M . 
The CHE does not occur any more for more massive stars 
(M ZAMS > 190 M Q ). 

4. The rotating models are kept rotating almost rigidly until the 
end of core helium burning, which results from magnetic 
torques due to the Spruit-Tayler dynamo (Sect. |3.3) . The 
surface can reach the critical rotation during the evolution, 
and undergo mass loss via a centrifugally driven wind, 
which may form a decretion disk, if the initial rotational 
velocity is sufficiently high (Fig. HJ. In this case, mass loss 
occurs mostly on the main sequence, in which the wind 
material is mainly composed of hydrogen and helium, if 
the star does not undergo the CHE. With CHE, however, 
significant mass loss can occur even during the post-main 
sequence phases, which may lead to significant enrichment 
of CNO elements in the surrounding medium. 

5. Chemical mixing by convection and/or rotation may lead to 
an increase of the neutron excess in the stellar core during 
the helium-burning phase (Sect. 4.3 1. The most notable 



enhancement is observed in the CHE models. However, this 
effect is still weak compared to the expected increase of the 
neutron excess during the more advanced nuclear burning 
stages, and it would have little impact in the nucleosynthesis 
ofPISNe. 

The mass limits for PISN may decrease with rota- 
tion, as the core size increases with rotationally in- 
duced mixing (Sect. |6.2| Fig. 12 1. In particular, stars of 
84 < Mzams/M < 190 may produce hydrogen-free PISNe 
(type Ibc), via CHE. Such PISNe may look like the type 
Ic SN 2007bi ( |Gal-Yam et al~||2009l ), which is believed to 
originate from the pair instability. 

If a Pop III star does not undergo CHE, the core cannot 
retain enough angular momentum to produce an explosion 
powered by rotation like a jet-driven supernova or gamma- 
ray burst (Sects. |3.3| and 6.3 1. With CHE, conditions for 



gamma-ray burst (GRB) progenitors may be fulfilled for 
13 < M Z ams/M ~ 84. All of our GRB progenitor models 
suffer mass loss caused by rotation until the end. This may 
lead to formation of a complicated aspherical circumstellar 
structure, which should significantly deviate from the typical 
density profile around a hot star (p oc r~ 2 ). On the other 
hand, super-collapsars from M > 250 M Q are not likely to 
occur. 

8. For the mass range of 56 < Mzams/Mq < 84 in the CHE 
regime, both pulsational PISN and GRB may occur one after 
another from the same progenitor (Fig. |T2]>. In this case, the 
GRB jet would penetrate a very dense and massive medium 
that has been formed shortly before core-collapse by the 
ejecta of the pulsational PISN, and an extremely bright af- 
terglow may be produced as a consequence. 
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Appendix A: Kippenhan diagram 
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Table 3. Total energy (E), total number (N), and average wavelength (A) of the ionizing photons for hydrogen (H), neutral helium 
(He) and singly ionized helium (He + ) during the evolution for each model sequence. 





E(th) 


N(th) 


My») 


E(y H e) 


N( 7He ) 




E(7He+) 


N(y H e+) 


^(THc+) 




[erg] 




[A] 


[erg] 




[A] 


[erg] 




[A] 


mlOvkOO 


1.69E+52 


5.51E+62 


648.51 


3.59E+51 


7.62E+61 


421.82 


1.22E+49 


1.30E+59 


211.04 


ml0vk02 


1.61E+52 


5.28E+62 


653.63 


3.23E+51 


6.89E+61 


423.67 


9.31E+48 


9.91E+58 


211.43 


ml0vk04 


1.36E+52 


4.55E+62 


667.55 


2.32E+51 


5.01E+61 


428.98 


3.87E+48 


4.14E+58 


212.65 


ml0vk06 


1.72E+52 


5.73E+62 


663.56 


3.08E+51 


6.64E+61 


427.52 


6.06E+48 


6.47E+58 


212.26 


ml0vk08 


1.56E+52 


5.24E+62 


667.37 


2.68E+51 


5.78E+61 


428.85 


4.59E+48 


4.91E+58 


212.54 


ml5vk00 


3.46E+52 


1.06E+63 


608.94 


1.05E+52 


2.15E+62 


405.57 


1.24E+50 


1.30E+60 


207.25 


ml5vk02 


3.61E+52 


1.11E+63 


612.28 


1.07E+52 


2.19E+62 


406.87 


1.16E+50 


1.21E+60 


207.59 


ml5vk04 


3.46E+52 


1.09E+63 


625.25 


9.20E+51 


1.91E+62 


412.56 


6.67E+49 


7.01E+59 


208.95 


ml5vk06 


1.10E+53 


2.88E+63 


518.98 


5.56E+52 


1.00E+63 


358.72 


5.81E+51 


5.49E+61 


187.98 


m20vk00 


5.66E+52 


1.67E+63 


586.19 


2.03E+52 


4.04E+62 


395.80 


4.20E+50 


4.33E+60 


204.93 


m20vk02 


5.70E+52 


1.69E+63 


588.77 


2.01E+52 


4.02E+62 


396.95 


3.91E+50 


4.03E+60 


205.25 


m20vk03 


7.69E+52 


2.25E+63 


580.70 


2.86E+52 


5.68E+62 


395.21 


5.86E+50 


6.05E+60 


205.22 


m20vk04 


1.63E+53 


4.06E+63 


495.57 


9.05E+52 


1.60E+63 


351.31 


1.09E+52 


1.02E+62 


186.83 


m20vk06 


1.55E+53 


3.89E+63 


499.29 


8.48E+52 


1.50E+63 


352.48 


9.95E+51 


9.37E+61 


187.08 


m30vk00 


1.05E+53 


2.95E+63 


558.54 


4.45E+52 


8.58E+62 


383.51 


1.63E+51 


1.66E+61 


202.00 


m30vk02 


1.06E+53 


3.02E+63 


563.56 


4.39E+52 


8.49E+62 


384.82 


1.53E+51 


1.56E+61 


202.31 


m30vk03 


2.61E+53 


6.15E+63 


468.83 


1.60E+53 


2.75E+63 


342.02 


2.25E+52 


2.11E+62 


186.35 


m30vk04 


2.56E+53 


6.11E+63 


473.67 


1.55E+53 


2.68E+63 


343.86 


2.10E+52 


1.97E+62 


186.46 


m30vk06 


2.42E+53 


5.81E+63 


476.28 


1.45E+53 


2.51E+63 


344.72 


1.95E+52 


1.83E+62 


186.42 


m60vk00 


2.58E+53 


6.76E+63 


521.37 


1.31E+53 


2.41E+63 


366.30 


9.02E+51 


8.98E+61 


197.72 


m60vk02 


4.09E+53 


1.04E+64 


503.97 


2.24E+53 


4.09E+63 


363.14 


1.63E+52 


1.63E+62 


197.76 


m60vk03 


5.41E+53 


1.20E+64 


442.55 


3.62E+53 


6.04E+63 


331.81 


6.03E+52 


5.63E+62 


185.67 


m60vk04 


5.27E+53 


1.18E+64 


446.09 


3.48E+53 


5.83E+63 


333.05 


5.70E+52 


5.32E+62 


185.62 


m60vk06 


4.93E+53 


1.11E+64 


448.25 


3.24E+53 


5.44E+63 


333.95 


5.22E+52 


4.88E+62 


185.70 


mlOOvkOO 


4.63E+53 


1.18E+64 


504.36 


2.52E+53 


4.54E+63 


357.71 


2.26E+52 


2.23E+62 


195.39 


ml00vk02 


5.20E+53 


1.23E+64 


469.73 


3.23E+53 


5.67E+63 


348.89 


3.44E+52 


3.37E+62 


194.88 


ml00vk03 


9.12E+53 


1.98E+64 


432.21 


6.29E+53 


1.04E+64 


327.59 


1.12E+53 


1.05E+63 


185.58 


ml00vk04 


8.79E+53 


1.92E+64 


434.24 


6.02E+53 


9.95E+63 


328.32 


1.06E+53 


9.94E+62 


185.56 


ml00vk06 


8.31E+53 


1.82E+64 


435.35 


5.67E+53 


9.38E+63 


328.70 


9.98E+52 


9.31E+62 


185.45 


ml50vk00 


7.03E+53 


1.75E+64 


494.17 


3.97E+53 


7.04E+63 


352.50 


4.11E+52 


4.01E+62 


193.98 


ml50vk02 


8.88E+53 


2.19E+64 


489.55 


5.11E+53 


9.08E+63 


353.13 


5.17E+52 


5.06E+62 


194.31 


ml50vk03 


1.30E+54 


2.87E+64 


440.71 


8.78E+53 


1.48E+64 


334.53 


1.32E+53 


1.26E+63 


190.80 


ml50vk04 


1.28E+54 


2.79E+64 


431.75 


8.87E+53 


1.46E+64 


327.87 


1.57E+53 


1.47E+63 


186.17 


ml50vk06 


1.23E+54 


2.66E+64 


430.85 


8.49E+53 


1.40E+64 


327.17 


1.52E+53 


1.42E+63 


185.76 


m200vk00 


9.35E+53 


2.30E+64 


489.40 


5.37E+53 


9.46E+63 


349.92 


5.96E+52 


5.80E+62 


193.22 


m200vk02 


1.41E+54 


3.45E+64 


485.24 


8.29E+53 


1.47E+64 


351.78 


8.54E+52 


8.35E+62 


194.33 


m200vk04 


1.59E+54 


3.65E+64 


454.98 


1.04E+54 


1.78E+64 


341.18 


1.33E+53 


1.29E+63 


193.04 


m200vk06 


1.59E+54 


3.53E+64 


442.45 


1.07E+54 


1.81E+64 


335.60 


1.56E+53 


1.50E+63 


191.39 


m250vk00 


1.11E+54 


2.70E+64 


484.51 


6.48E+53 


1.13E+64 


347.65 


7.59E+52 


7.36E+62 


192.68 


m250vk02 


1.36E+54 


3.31E+64 


483.27 


8.03E+53 


1.41E+64 


349.36 


8.87E+52 


8.64E+62 


193.59 


m250vk04 


1.89E+54 


4.40E+64 


462.38 


1.20E+54 


2.08E+64 


344.63 


1.43E+53 


1.39E+63 


193.83 


m250vk06 


1.76E+54 


4.10E+64 


462.31 


1.12E+54 


1.94E+64 


344.48 


1.34E+53 


1.30E+63 


193.73 


m300vk00 


1.35E+54 


3.30E+64 


484.31 


7.91E+53 


1.38E+64 


347.22 


9.41E+52 


9.11E+62 


192.43 


m300vk02 


1.89E+54 


4.61E+64 


485.29 


1.11E+54 


1.96E+64 


351.77 


1.15E+53 


1.12E+63 


194.04 


m300vk04 


2.00E+54 


4.69E+64 


465.91 


1.26E+54 


2.19E+64 


346.47 


1.43E+53 


1.40E+63 


194.14 


m500vk00 


2.15E+54 


5.20E+64 


480.96 


1.27E+54 


2.21E+64 


345.26 


1.59E+53 


1.53E+63 


191.81 


m500vk04 


2.45E+54 


6.00E+64 


487.16 


1.42E+54 


2.52E+64 


351.39 


1.49E+53 


1.45E+63 


194.33 
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Table 4. Yields of stellar winds from rotating models, in units of solar mass 





'H 


4 He 


12 C 


13 C 


14 N 


16q 


18q 


22 Ne 


ml0vk02 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


ml0vk04 


1.11E-01 


4.47E-02 


2.56E-22 


8.29E-23 


3.79E-20 


1.11E-21 


6.85E-28 


8.26E-58 


ml0vk04 


4.05E-01 


1.83E-01 


9.18E-17 


2.97E-17 


1.29E-14 


3.50E-16 


1.88E-22 


1.02E-22 


ml0vk06 


6.37E-01 


2.55E-01 


1.24E-20 


4.02E-21 


1.82E-18 


5.26E-20 


3.20E-26 


6.37E-36 


ml5vk02 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


ml5vk04 


3.79E-01 


1.67E-01 


7.36E-14 


2.37E-14 


9.36E-12 


2.39E-13 


1.14E-19 


2.53E-20 


ml5vk06 


5.32E-01 


2.45E+00 


2.46E-02 


2.08E-03 


9.79E-03 


4.76E-01 


3.59E-06 


2.47E-04 


m20vk02 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


m20vk03 


3.58E-01 


2.31E-01 


2.06E-12 


6.64E-13 


2.29E-10 


5.21E-12 


2.00E-18 


6.26E-18 


m20vk04 


3.80E-02 


2.03E+00 


6.08E-02 


1.84E-03 


7.22E-03 


1.36E-01 


1.98E-06 


2.35E-04 


m20vk06 


6.90E-01 


3.23E+00 


6.28E-02 


5.29E-03 


1.19E-02 


8.07E-01 


3.80E-05 


8.70E-04 


m30vk02 


1.77E-01 


6.43E-02 


3.34E-13 


1.08E-13 


3.63E-11 


8.40E-13 


3.24E-19 


6.72E-19 


m30vk03 


5.61E-02 


3.92E+00 


6.22E-02 


7.92E-03 


2.52E-02 


1.47E+00 


2.30E-05 


2.28E-03 


m30vk04 


1.41E-01 


4.90E+00 


4.94E-02 


7.68E-03 


2.24E-02 


1.43E+00 


1.30E-05 


1.88E-03 


m30vk06 


1.26E+00 


4.98E+00 


1.92E-02 


6.81E-04 


2.06E-02 


1.40E+00 


1.53E-06 


2.23E-03 


m60vk02 


1.52E+00 


2.06E+00 


1.61E-10 


5.17E-11 


1.43E-08 


2.83E-10 


8.85E-17 


9.24E-16 


m60vk03 


2.08E-01 


9.51E+00 


9.72E-02 


1.01E-02 


3.45E-02 


3.34E+00 


3.59E-06 


4.36E-03 


m60vk04 


1.02E+00 


1.08E+01 


9.62E-02 


8.32E-03 


2.21E-02 


3.26E+00 


3.11E-06 


3.64E-03 


m60vk06 


3.87E+00 


1.72E+01 


1.36E-01 


1.80E-02 


4.20E-02 


4.39E+00 


9.97E-06 


5.64E-03 


ml00vk02 


3.73E+00 


3.75E+00 


3.12E-10 


9.96E-11 


2.67E-08 


4.49E-10 


1.49E-16 


1.77E-15 


ml00vk03 


9.12E-01 


1.74E+01 


6.18E-01 


1.53E-02 


3.30E-02 


4.53E+00 


2.54E-05 


2.90E-03 


ml00vk04 


2.97E+00 


1.84E+01 


3.20E-01 


1.96E-02 


4.04E-02 


3.82E+00 


1.04E-05 


2.78E-03 


ml00vk06 


6.80E+00 


1.87E+01 


3.27E-01 


2.72E-03 


4.91E-02 


3.71E+00 


3.97E-07 


2.49E-03 


ml50vk02 


6.62E+00 


5.88E+00 


5.57E-10 


1.77E-10 


4.73E-08 


7.85E-10 


2.60E-16 


3.64E-15 


ml50vk03 


4.65E+00 


2.34E+01 


1.61E-08 


5.14E-09 


1.11E-06 


1.34E-08 


4.44E-15 


5.52E-11 


ml50vk04 


7.24E+00 


2.67E+01 


1.35E+00 


1.56E-02 


4.88E-02 


2.05E+00 


4.60E-07 


3.25E-04 


ml50vk06 


1.22E+01 


2.76E+01 


1.81E+00 


1.71E-02 


4.66E-02 


3.23E+00 


8.78E-07 


6.80E-04 


m200vk02 


7.22E+00 


8.82E+00 


1.25E-09 


4.02E-10 


1.01E-07 


1.72E-09 


4.98E-16 


9.36E-15 


m200vk04 


1.38E+01 


2.62E+01 


5.80E-09 


1.86E-09 


4.54E-07 


5.49E-09 


1.67E-15 


4.82E-14 


m200vk06 


1.96E+01 


3.39E+01 


1.02E-08 


3.27E-09 


7.48E-07 


8.42E-09 


2.78E-15 


8.95E-14 


ml50vk02 


1.16E+01 


9.13E+00 


1.06E-09 


3.35E-10 


8.26E-08 


1.27E-09 


3.98E-16 


6.60E-15 


m250vk04 


1.99E+01 


3.08E+01 


6.70E-09 


2.17E-09 


5.19E-07 


4.11E-09 


2.60E-15 


5.75E-14 


m250vk06 


2.75E+01 


3.24E+01 


5.65E-09 


1.81E-09 


4.59E-07 


3.77E-09 


2.27E-15 


4.82E-14 


m300vk02 


1.25E+01 


1.13E+01 


1.42E-09 


4.56E-10 


1.21E-07 


2.09E-09 


6.06E-16 


1.16E-14 


m300vk04 


2.80E+01 


3.63E+01 


6.13E-09 


1.98E-09 


4.98E-07 


4.21E-09 


2.49E-15 


5.46E-14 


m500vk04 


5.53E+01 


3.14E+01 


1.06E-08 


1.18E-09 


2.64E-07 


2.58E-08 


1.36E-15 


3.01E-13 
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Table 5. Yields in the last calculated model for each model sequence, in units of solar mass. The wind yields are not included here. 
?7 C denotes the neutron excess at the center of the last model. X max ( 14 N)is the maximum mass fraction of 14 N that has been achieved 
during the evolution. 



"TS 



H 



4 He 



N 



o 



O 



Ne 



X 



4.90E+00 
4.94E+00 
5.04E+00 
4.17E+00 
4.12E+00 



2.92E+00 
2.95E+00 
2.96E+00 
3.19E+00 
3.11E+00 



1.84E-01 
1.80E-01 
9.21E-02 
1.74E-01 
1.50E-01 



2.44E-09 
2.39E-09 
1.39E-09 
2.10E-09 
1.79E-09 



6.13E-08 
5.35E-08 
2.83E-08 
4.49E-08 
3.59E-08 



1.04E+00 
1.05E+00 
8.84E-01 
1.06E+00 
9.42E-01 



8.18E-09 
8.18E-09 
7.27E-09 
9.97E-09 
7.18E-09 



1.53E-07 
1.96E-07 
8.48E-08 
1.95E-07 
1.30E-07 



4.69E-08 
4.45E-08 
3.95E-08 
4.24E-08 
3.16E-08 



6.98E-07 
6.47E-07 
6.26E-07 
6.25E-07 
5.86E-07 



6.75E+00 
6.46E+00 
6.03E+00 
1.99E-20 
8.26E+00 
8.08E+00 
5.67E+00 
1.97E-20 
1.05E-20 



4.35E+00 
4.49E+00 
4.70E+00 
6.24E-01 
5.55E+00 
5.96E+00 
6.69E+00 
4.55E-01 
5.61E-01 



3.57E-01 
3.21E-01 
3.32E-01 
1.99E-01 
1.04E+00 
1.37E+00 
1.80E+00 
9.89E-01 
2.51E-01 



2.82E-09 
2.63E-09 
5.16E-09 
1.55E-04 
4.57E-07 
4.85E-05 
7.00E-07 
2.80E-04 
3.80E-04 



1.08E-07 
8.32E-08 
1.05E-07 
6.13E-04 
5.69E-08 
8.16E-04 
1.21E-07 
1.30E-03 
1.24E-03 



2.12E+00 
1.71E+00 
2.13E+00 
9.79E+00 
4.69E+00 
4.16E+00 
5.12E+00 
1.33E+01 
1.29E+01 



7.02E-08 
7.69E-08 
9.26E-08 
6.09E-06 
1.36E-07 
1.47E-06 
2.51E-07 
5.44E-04 
4.66E-04 



3.23E-07 
3.81E-07 
4.62E-07 
4.81E-03 
2.72E-06 
3.17E-06 
4.30E-06 
3.67E-03 
2.86E-03 



6.56E-08 
8.33E-08 
7.95E-08 
6.90E-05 
1.36E-07 
1.22E-07 
1.41E-07 
1.58E-04 
1.25E-04 



1.10E-06 
8.50E-07 
8.42E-07 
1.43E-02 
9.60E-07 
1.06E-03 
1.21E-06 
1.40E-02 
1.16E-02 



1.07E+01 
1.00E+01 
1.85E-20 
2.03E-20 
1.96E-20 



7.96E+00 
8.00E+00 
4.18E-01 
4.98E-01 
2.76E-01 



3.37E-01 
5.22E-01 
3.99E-01 
3.45E-01 
4.29E-01 



1.13E-05 
2.90E-08 
2.96E-04 
3.15E-04 
2.60E-05 



2.58E-04 
2.88E-07 
1.27E-03 
1.17E-03 
7.65E-04 



9.21E+00 
1.04E+01 
2.08E+01 
1.98E+01 
1.91E+01 



1.44E-09 
1.17E-07 
4.92E-04 
4.78E-04 
1.99E-04 



3.49E-06 
6.13E-06 
2.37E-03 
2.40E-03 
1.50E-03 



1.48E-06 
1.10E-06 
1.69E-04 
1.41E-04 
1.70E-04 



1.35E-02 
1.57E-04 
1.48E-02 
1.61E-02 
1.22E-02 



1.67E+01 
5.40E+00 
2.18E-20 
2.68E-20 
2.67E-20 



1.41E+01 
1.44E+01 
2.81E-01 
3.02E-01 
2.90E-01 



1.75E+00 
2.05E+00 
9.48E-01 
9.25E-01 
7.89E-01 



7.61E-08 
3.34E-05 
6.44E-05 
5.45E-05 
6.66E-05 



6.85E-07 
9.21E-05 
2.31E-04 
1.41E-04 
1.43E-04 



2.49E+01 
3.09E+01 
3.84E+01 
3.70E+01 
3.49E+01 



8.30E-07 
1.59E-05 
3.00E-05 
1.94E-05 
2.60E-05 



5.29E-06 
5.34E-06 
6.88E-04 
6.07E-04 
7.60E-04 



1.86E-06 
8.26E-07 
1.16E-04 
9.32E-05 
1.05E-04 



2.38E-04 
1.47E-04 
1.11E-02 
7.88E-03 
7.91E-03 



2.52E+01 
4.25E+00 
2.01E-20 
2.00E-18 
4.52E-19 



2.29E+01 
2.17E+01 
4.51E-01 
5.97E-01 
5.58E-01 



2.68E+00 
3.84E+00 
2.71E+00 
2.26E+00 
2.50E+00 



1.77E-06 
6.25E-07 
8.30E-05 
2.54E-04 
5.10E-05 



2.87E-06 
3.29E-06 
1.64E-04 
5.43E-04 
7.45E-04 



4.32E+01 
5.42E+01 
6.09E+01 
5.91E+01 
5.65E+01 



3.50E-07 
5.16E-06 
3.37E-06 
3.50E-05 
4.77E-05 



7.31E-06 
9.58E-05 
3.84E-04 
6.49E-04 
5.43E-04 



6.14E-06 
7.13E-07 
4.80E-05 
5.57E-05 
5.17E-05 



9.96E-03 
1.75E-04 
7.42E-03 
9.11E-03 
1.37E-02 



3.24E+01 
9.37E+00 
8.32E-02 
1.56E-18 
1.28E-18 
3.97E+01 
1.42E+01 
1.96E+00 
1.75E-01 



3.54E+01 
2.72E+01 
9.86E+00 
7.19E-01 
3.99E-01 
5.30E+01 
4.48E+01 
2.56E+01 
1.29E+01 



5.27E+00 
5.12E+00 
8.78E+00 
7.96E+00 
6.83E+00 
8.49E+00 
7.04E+00 
7.92E+00 
9.54E+00 



1.52E-06 
7.97E-07 
1.56E-04 
3.48E-04 
9.79E-05 
1.51E-03 
2.89E-06 
1.96E-06 
3.91E-05 



3.49E-05 
5.89E-05 
1.15E-03 
1.45E-03 
3.37E-04 
3.71E-04 
6.54E-05 
5.44E-06 
5.10E-04 



6.67E+01 
7.96E+01 
8.82E+01 
8.86E+01 
8.36E+01 
8.54E+01 
9.74E+01 
1.02E+02 
1.03E+02 



1.18E-08 
3.13E-05 
4.89E-04 
2.23E-06 
1.11E-06 
2.58E-11 
1.90E-05 
6.20E-05 
6.00E-04 



1.26E-05 
2.44E-04 
5.23E-05 
1.53E-04 
9.69E-05 
1.83E-05 
7.52E-05 
1.60E-04 
5.11E-05 



2.23E-06 
1.19E-05 
1.35E-06 
1.20E-05 
1.72E-05 
6.36E-07 
1.67E-06 
4.76E-07 
1.02E-06 



1.10E-03 
3.72E-03 
1.62E-03 
1.40E-02 
1.07E-02 
3.44E-03 
4.46E-04 
7.90E-05 
3.96E-04 



4.27E+01 
2.19E+01 
1.28E+00 
2.06E+00 



6.66E+01 
5.31E+01 
2.67E+01 
2.90E+01 



1.06E+01 
7.81E+00 
1.02E+01 
9.13E+00 



1.24E-03 
7.93E-06 
7.17E-05 
9.13E-06 



7.68E-05 
6.48E-05 
4.15E-04 
1.94E-04 



1.10E+02 
1.18E+02 
1.31E+02 
1.22E+02 



3.26E-10 
1.12E-05 
3.92E-09 
1.48E-04 



2.37E-05 
2.94E-04 
3.84E-06 
1.11E-05 



7.10E-07 
3.96E-06 
6.69E-07 
6.83E-07 



1.70E-04 
1.76E-03 
2.84E-04 
3.52E-04 



5.33E+01 
3.00E+01 
4.19E+00 



8.25E+01 
8.01E+01 
4.20E+01 



1.27E+01 
9.45E+00 
9.89E+00 



1.50E-02 
2.62E-05 
2.21E-06 



1.27E-03 
1.52E-04 
1.24E-04 



1.26E+02 
1.27E+02 
1.43E+02 



4.67E-10 
8.47E-05 
1.74E-04 



4.29E-05 
3.19E-05 
9.39E-06 



6.03E-07 
4.10E-07 
5.29E-07 



9.36E-03 
1.40E-04 
2.18E-04 



8.46E+01 
2.63E+01 



1.45E+02 
9.08E+01 



2.22E+01 
1.47E+01 



3.41E-01 
5.33E-04 



1.17E-01 
3.43E-04 



1.97E+02 
2.17E+02 



2.30E-09 
6.56E-11 



3.57E-05 
4.34E-06 



1.95E-07 
3.92E-06 



8.59E-03 
1.78E-03 
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